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Abstract 

Generalisations of the 't Hooft-Polyakov monopole which can exhibit re- 
pulsion only, attraction only, and both attraction and repulsion, between like 
monopoles, are studied numerically. The models supporting these solitons are 
SO (3) gauged Higgs models featuring Skyrme-like terms. 
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1 Introduction 



It is well known that 't Hooft-Polyakov 0] monopoles of like magnetic charge 
repel. For all finite values of the Higgs mass, or the coupling constant of the Higgs 
self-interaction potential, like monopoles exhibit a single repulsive phase only. This 
is understood as follows: In the Prasad-Sommerfield (PS) limit, the Higgs field is 
massless and hence mediates a long range attractive force that cancels the long range 
repulsive magnetic force of the U(l) field, exactly. In this limit, all components of 
the stress tensor vanish j|] resulting in the vanishing of inter-monopole forces. The 
force between like monopoles in the BPS limit was studied in detail by Manton || 
and Nahm ||, who showed that it decreases faster than any inverse power. In the 
presence of the usual Higgs potential however, the Higgs field becomes massive and as 
a result decays exponentially. Consequently the long range magnetic field dominates 
at large distances, leading to the repulsion of like monopoles of the Georgi-Glashow 
(GG) model. This was also concluded by Goldberg et al. f7|, using the time-rate of 
change of the stress tensor for the field configuration of two exponetially localized 
monopoles, situated apart at a distance much larger than the sizes of the monopole 
cores. Recently Kiselev and Shnir |5| also have shown this employing the fluctuations 
and zero modes around the monopoles. That there is only one repulsive phase has 
been verified in , by calculating the energy per unit topological charge for a charge- 
2 and charge-1 solution numerically for all values of the Higgs potential coupling 
constant. The procedure in Ref. || is exactly analogous to that adopted by Jacobs 
and Rebbi |TU] in their investigation of attractive and repulsive phases in the Abelian 
Higgs model. 

Implementing the procedure of Ref. |TI| in our 3 dimensional context necessitates 
the numerical solution of the partial differential equations of the axially symmetric 
fields. This is because the charge-2 monopole in 3 dimensions is not spherically sym- 



metric, unlike in the case of the 2 dimensional context of the Abelian Higgs model [10 



and of the hierarchy of Abelian Higgs models [H], where vorticity-n solutions are ra- 
dially symmetric and hence result from ordinary differential equations. Thus, one of 
the most important technical ingredients of this work is the imposition of the axial 
symmetry on the 50(3) gauge field and the algebra valued Higgs, for which we use 



the Ansatz of Rebbi and Rossi [12] 



The purpose of the present work is to show that certain generalisations of the 
monopoles of the GG model, exhibit (a) a repulsive phase only like the usual monopole, 
(b) an attractive only, and (c) both attractive and repulsive phases. This is done by 
integrating the Euler-Lagrange equations arising from the static Hamiltonians of the 
new models, numerically, in exact analogy with the procedures used in Refs. fl(]. U . 
The models under consideration, just like the GG model, are 3 dimensional £0(3) 
Higgs models. These models are constructed from the hierarchy of 3 dimensional 
£0(3) Higgs models, obtained by dimensional descent from the 4p dimensional hier- 
archy of chiral SO(Ap) Yang-Mills models [13] Q. Here, we will restrict to the p — 1 



1 There are similar hierarchies of d dimensional SO(d) Higgs models with d < Ap Il4 
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and the p = 2 members of this hierarchy, of which the p = 1 member is nothing else 
than the PS limit of the GG model. Thus the main new ingredient in the models 
presented here will be the p = 2 model, whose spherically symmetric monopol solu- 
tions were studied numerically in our previous work to which we shall refer for 
certain details later. 



In Ref. ||15|| , we have discussed the charge-1 solutions of the p = 2 model in some 
detail, especially with regard to its properties relative to the GG and the Abelian 
Higgs models. There we have pointed out that the p = 2 model does not support self- 
dual solutions. However, for certain values of the coupling constants it nevertheless 
supports a solution whose energy is quantitatively extremely close to the Bogomol'nyi 
lower bound, and correspondingly, an estimate of some components of the stress ten- 
sor indicates that the force on the soliton is also extremely small. We described such 
a solution as being "almost self-dual" and we designated the corresponding values 
of the coupling constants as critical. This is in accordance with the usual nomen- 
clature whereby the Bogomol'nyi inequalities are saturated for the critical values of 



the coupling constants in a given model. It is important to recall |T4], [15| that the 
critical dimensionless coupling constants of the p = 2 model take nonzero values. In 
this sense the p = 2 model is similar to the Abelian Higgs model, which supports a 
self-dual solution for a finite value of the dimensionless coupling constant, and it is in 
contrast to the GG model for which the critical value of the (Higgs) coupling constant 
vanishes. Hence it was argued in Ref. [15| that like the Abelian Higgs model, the 3 
dimensional p = 2 SO (3) Higgs model, supports both attractive and repulsive phases. 

To put all this in perpective, we note that there is a difference between the p = 2 
model on the one hand and the p = 1 and the abelian Higgs model on the other hand. 
In the latter, the Bogolmol'nyi inequality is saturated at the critical coupling con- 
stants for all topological charges n. Thus the critical values of the coupling constants 
characterise a property of the model itself. For the p = 2 model we only know that 
for the charge-1 solution the energy is close to the Bogolmol'nyi bound at the critical 
coupling constants. If these coupling constants characterise the model, and not only 
a single solution, then the energy of all charge-n solutions should be near their Bo- 
golmol'nyi bounds at the same values of the coupling constants. So we do not specify 
quantatively what we mean by 'near the Bogolmol'nyi bound', there might be a sub- 
set of values of the coupling constants which can be called critical. For the charge-2 
solution of the p = 2 model we expect that the energy will be near the Bogolmol'nyi 
bound for coupling constants close to the critical coupling constants determined by 
the investigation of the charge-1 solutions. If this is the case, attractive and repulsive 
phases possibly exist with a crossover at values of the coupling constants near the 
critical values. 

It is one of the aims of the present work to verify this hypothesis. This will be 
done by straighforwardly integrating the Euler-Lagrange equations numerically on the 
charge-1 and charge-2 sectors. Thus the argument based on the almost self- duality of 
the p = 2 monopoles, that makes the existence of the two phases plausible, will not 



be further elaborated here, and we refer to Refs. 11, 15 for that. 
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In general, our aim is to present all the generalised 50(2) Higgs models with 
Hamiltonian densities carrying dimensions up to L~ 8 , which includes of course the 
GG model with dimension L~ 4 . Note that the limitation we impose is a consequence of 
limiting our study to the p = 2 model, whose Hamiltonian density has the dimensions 
L~ 8 . The inclusion of the GG or the p — 1 model results in the presence of the usual 
Yang-Mills term and hence in the possibility of describing free gauge fields, absent 
in the pure p = 2 (and p > 2) model(s). We note that the p > 2 models are purely 
Skyrme-like gauged Higgs systems [0] in the sense that they consist of 2pth powers 
of the curvature and covariant derivative, suitably antisymmetrised so that only the 
squares of velocity-fields appear. In this sense, the composite models consisting of 
the p = 1 model augmented by p > 2 models, are quite analogous to the Skyrme 
model [T7J where the latter terms play the role of the Skyrme term of the 0(4) sigma 
model. 

It is worthwhile making some clarifications concerning the construction of the 



composite models. We know from the results and the arguments of Refs. [|TT], [15 
that we would expect a phase change from attractive to repulsive, to take place at 
the almost self-dual configurations of the systems. The composite models however 
are even less self-dual than the p — 1 and the p = 2 modes on their own, since 
the Bogolmol'nyi equations are even more overdetermined then those of the p = 2 
model on its own and hence the solutions of the composite models are more remote 
from saturating the Bogolmol'nyi bound. (Clearly, exact self-dual solutions can be 
achieved if we suppress the entire p = 2 content of the model, but then there will be no 
interactions between the monopoles. Since our aim is to produce a (composite) model 
accomodating interactions, we will resist on retaining some or all of the components 
of the p = 2 model.) We will investigate three possibilities, (a) when the p = 1 system 
dominates, which can be regarded as a p = 1 model extended with a p = 2 Skyrme- 
term, (b) when the two constituents of the composite model are equally important, 
and (c) when the p = 2 model dominates, which yields the pure p = 2 model whose 
spherically symmetric solutions were studied in Ref. |]15 |. We shall study all these 
three cases in this work. 

In Section 2 we define the static Hamiltonians of the various models we will study 
and proceed in Section 3 to state the axially symmetric Ansatz of Rebbi and Rossi [fHql , 
leading to the residual 2 dimensional subsytems of the models introduced in the 
previous section. In Sections 4, 5 and 6, we present the results of the numerical 
computations, respectively for the extended p = 1 models where the latter dominates, 
the composite p = 1 and p = 2 system, and the pure p = 2 system which is the limiting 
case of the composite system where the p = 2 dominates. Our choice for this order 
of presentation is the result of the successively increasing complication involved in 
the numerical computations involved in each of these cases, so that we have chosen 
to tackle the most easily accessible cases first. Section 7, finally, is devoted to the 
summary and discussion of the results. 
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2 The Models 



Since various aspects of the 3 dimensional hierarchy of SO (3) gauged Higgs models 



are discussed in detail in Refs. 14, U5L [L6|, we shall restrict ourselves here to the 



definitions of the Hamiltonian densities only. We refer to Refs. p|, ITB] for their 



construction, as well as the topological inequalities, via dimensional descent. For 
the detailed normalisation of the topological charge and the asymptotic properties 
required by the finite energy condition and analyticity, we refer to Ref. Jl5| . 
Below, we shall denote the gauge and Higgs fields 

A = ~\^>a , <P = -]fVa (1) 

in an antihermitian notation. The resulting antihermitian curvature and covariant 
derivative are then 

Fij = diAj - OjAi + [A l} Aj] , (2) 
A0 = d^ + [A,# (3) 

2.1 The p = 1 model 

This is the PS limit of the GG model, characterised by the static Hamiltonian 

= Tr(l /••,,/••'' + ^DrfDty . (4) 

The reason that we have distinguished between covariant and contravariant indices is 
that in what follows we shall have occasion to use a spherical (not Cartesian) basis. 



The system supports self-dual solutions saturating the topological lower bound |18 
which equals the magnetic charge. This model is extensively studied since a long time 
and will not be considered in isolation here. 

2.2 The p = 2 model 



The static Hamiltonian of this model is 15 



where each of the four densities H.^ are defined by 

Hf = Tr^D^HF^ , (6) 

Hf ] = -6Ti({S,F 3 } + [D t( p,D 3( p})({S,F^} + [D^,D^]) , (7) 

Uf = -27Tr[{S,A0}{S,£>V}] , (8) 

Hf ] = 54TrS 4 . (9) 
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In we have used [ijk] to denote cyclic symmetry in the indices k, {A, B} 

denotes an anticommutator, and 



S = -( V 2 + r). (10) 

The reason that we have distinguished between covariant and contravariant indices 
in ©-(H) is the same as in (f|). 

In (H) we can set A 4 = 1 by scaling with an overall constant. (A 4 is reserved for 
use only when we construct a composite model by combining the p = 1 and p = 2 
models.) When the dimensionless coupling constants Ai,A 2 and A 3 are restricted to 
their critical values Ai = A2 = A3 = 1, then the topological inequalitiy 

J n {2 \\i = A 2 = A 3 = l)d 3 r > 47mr/ 5 , (11) 

where n is the topological charge, can be saturated by the corresponding Bogomol'nyi 
equations. The latter being overdetermined [|L4| they do not support nontrivial so- 
lutions. Nevertheless there exist nontrivial solutions, that do not satisfy the Bogo- 
mol'nyi equations, with energies lying above but numerically very close to the Bo- 



gomol'nyi lower bound. The spherically symmetric solution found in Ref. |15| is of 
this type, described as being almost self-dual. When all the coupling constants are 
greater than 1, i.e. A^ > 1, then the topological inequality (O) is preserved. When 
however any of the coupling constants takes a lower value than 1, i.e. A^ < 1, then 
the inequality ( |TT]) is replaced by 

J H {2) d 3 r > 47iX min n V 5 , (12) 

where \ m in is the smallest of A^, A = 1, 2, 3. 

^From the results in Ref. [|ll|], where we have seen that there is a crossover point 
between attractive and repulsive phases in the vicinity of the critical values of the 
dimensionless coupling constants of a system supporting almost self-dual solutions, 
we expect this to be the case also for the p = 2 model (Bf), which we know from the 



results of Ref. [15] that it does support almost self-dual solutions. 



2.3 The composite model 

By composite models we mean the system characterised by the sum of the static 
Hamiltonians H,( p \ always including the p — 1 member of the hierarchy. Here we 
are limiting ourselves to p — 1 and p — 2. In 2.3.1, we define the models we shall 
subsequently study numerically. In 2.3.2, we estimate the dependence of the energy 
on the coupling constants Xa in the region A^ 1. 

2.3.1 The generic case 

The Hamiltonian density of this model is 

H (c) =K (1) +f 4 K (2) , (13) 
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whose Bogomol'nyi lower bound is twice Aimr] if Aa > 1- 

In ([Uf), the dimensionless coupling constants Aa for all four A = 1,2, 3, 4, can take 
on any positive definite values, and even when any or all Aa vanish, the Bogomol'nyi 
lower bound does not sink to zero as for the p = 2 model in (|12"D, since the energy 
of the p = 1 model is still bounded by the usual magnetic charge. The fact that it 



is possible to set all but one of the Aa equal to zero in ([13]) leads us to describe (p"3|) 
as an extended p — 1 model. For example if we set all but one Aa in (|i~3|) equal to 
zero we have four special composite models, and if we set any two of them equal to 
zero we have six special composite models, etc. The first case is quite special. These 
four models are similar to the GG model in sofar as they describe a p = 1 model with 
only one additional term. This suggests that like the GG model, these four composite 
models support only one phase, i. e. an attractive only or a repulsive one only. We 
shall give a strongly suggestive analytic argument for this in 2.3.2 below. 

Since the Bogolmol'nyi bounds of the solutions to (|13"D are more substantially 
exceeded than in the p = 2 model, here we cannot invoke the presence of almost 
self-dual solutions. Thus in the generic case, where Aa are not considered to be small 
we cannot predict anything about the interaction ot the monopoles. The situation 
is different however, when all Aa <C 1. In these cases it is possible to examine this 
question systematically as we shall do in the following subsections. 

In order to determine whether there is a repulsive phase, an attractive phase or 



a crossover we follow Jacobs and Rebbi fllQl , and define the difference of the energy 



per unit topolgical charge for the charge-2 and charge- 1 solutions at fixed coupling 
constants A a, 

AE(\ A ) d M-^M- E V(\ A ), (14 ) 

where E^^Xa) is the energy of the charge-n solution. The monopoles repel when 
AE(X A ) > and attract when AE(X A ) < 0. If the quantity AE(X A ) vanishes, we 
find a crossover at Aa- 

2.3.2 Perturbed p = 1 models 

What is very interesting about the extended models is, that the behaviour of the 
solutions on the constants Aa can be analytically estimated for small values of the 
coupling constants, Aa <C 1- We shall restrict, at first, to the case where there is only 
one nonzero A a in US 2 '. 

Consider the energy integral 

E {n \X A ) = J (W« + J^V^aH^x, (15) 

as a function of Aa- The energy density functional in (|15D is the special case of the 
composite static Hamiltonian (O) in which all but a given coupling constant Aa are 



set equal to zero, and is evatuated for the topological charge-n solution described by 
the fields (4 (Aa) ,0( a ^). 
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We now insist that A a <C 1, and make a Taylor expansion of E^(\a 



E^ n \X A ) = £W(0) + A, 



d\, 



+ o(Ai). (16) 

X A =0 



Because Aa <C 1, we can describe the solution fields in the energy density functional 
in ( [TED as small perturbations of the known charge-ra PS || solution fields , < 



4 (Aa) = 4 (0) + Aa«„ (A ' 4) = {o) + \a V (17) 

in terms of the functions (dj, which in principle can be found. 

The first term in the Taylor expansion ([TBI) is recognised as the energy of the PS 
solution, and equals Aimrj. To calculate the second term, we differentiate the energy 



integral ([Tq) , which involves the differentiation of the functional in the integrand 
both with respect to its implicit dependence on A a through the fields A[ ,(j)^ A \ 
and its explicit dependence. The terms arising from the differentiation of the implicit 
dependence on Aa amount to subjecting the functional + 18 ^ 32 77~ 4 Aa7-^ 2 ) to the 
variational principle, which vanish by virtue of the Euler-Lagrange equations satisfied 
by the solution fields (|T7j). The result of the differentiation with respect to the explicit 



dependence on Aa then is the only contribution 



dE^ 5 



d\ A 18 x 32 
which evaluated at Aa = is 
dE^ 



dX. 



\ n 18 x 32 

Aa=0 



v- 4 n^[A^\^]d\ (it 



V~ 4 Jnf[Af\^]d^ (19) 



an integral evaluated strictly in terms of the PS solution fields (Af\ 0^). This is the 
sense in which the present calculation is described as being perturbative. 

Refining our notation for the PS solution fields (Af*\ <j)^) by labeling them further 
with the topological charge n as (Af'\n),<p^(n)), we substitute ( |i~9"l) in ([b 



AE = c A X A (20) 

5 l v 'H I ^ 2) [4 0) (2),^)(2)]d 3 x - / nf[M%U i0 Wxyi) 



CA 



18 x 32 



Clearly, when the quantity c A in fl20"D is positive (negative) the monopoles repel (at- 
tract). Moreover, if we assume that the dependence of the energies on Aa is monotonic, 
then this will be the only phase, like in the GG model. It is important however to 
stress that the above arguments do not apply to the GG model, namely, that the GG 
model cannot be treated as a perturbation of the p — 1 (BPS) system JJ. 

for the charge- 1 and the charge-2 solutions 



In table 1 the values of ] " ; 



n cL\a 



Aa=0 



and ca are given. Note, c A , A = 2,3,4 are of the same order of magnitude, whereas 
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ci is more than one order of magnitude larger. c% and c 2 are positive, leading to 
repulsion of like monopoles for the models with Ai 7^ or A2 7^ 0. For the models 
with A3 7^ or A3 7^ the constants C3 and C4 are negative. Hence these models have 
only attractive phases. 

The above perturbative argument can be systematically extended to the case of 
the full composite model with all or some of the A^ 7^ 0. The only difference there 
will be that the Taylor expansion will have more than one contribution (one from 
each nonzero Xa), 

4 

AE = "£c A X A . 

A=l 

(22) 

Indeed, it would be expected that there should be both phases, like in the generic 
case. If at least two of the constants ca have opposite signs then the model has 
attractive and repulsive phases. The change of the phases occurs for the values of the 
parameters A^ for which AE becomes zero. This is verified by our numerical studies 
for small values of Xa, reported later. 



3 Axial Symmetry 

Imposition of axial symmetry proceeds readily by the application of the Rebbi-Rossi 
Ansatz When subjected to axial symmetry, the 12 fields consisting of the 3 

components of the isospin-1 gauge field and the isospin-1 Higgs field are parametrised 
by 6 functions. It is instructive to describe [|T^, [Tj| these 6 funcions as a 2 component 



real field A^, with \i = 1, 2 labeling the 2 dimensional coordinate £ M = r,z, r = 
\/x 2 + y' 2 , and two complex scalar fields if and x- 111 this parametrisation, all Lorentz 
and SU(2) gauge invariant Hamiltonian densities H.^ reduce to an Abelian Higgs 
model in which A^ plays the role of a U(l) connection, and the two complex scalar 
fields ip and x pl a Y the roles of Higgs fields. The asymptotic behaviours required of 
these functions ip and x by the requirements of finite energy and analyticity are typical 
of Higgs fields, hence justifying the description of the axially symmetric 2 dimensional 
reduced subsytems of thep-th 5*0(3) gauged Higgs model in 3 dimensions, as residual 
Abelian two-Higgs models in 2 dimensions. (In particular, one of the complex fields, 
say ip actually is endowed with a Higgs type self-intereaction potential, and the two 
Higgs fields exhibit a Yukawa-type mutual interaction term.) Given that the residual 
subsytem exhibits a U(l) gauge symmetry, it is necessary to remove one of the 6 
residual fields by fixing this gauge. The natural choice we will make here, which has 
been used previously in Refs. PO, |2"Tf, is 



= °> ( 23 ) 
which we call a Coulomb gauge as before E0, |2T| . 
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In the following, we shall not describe the axially symmetric subsystem in terms 
of the functions (A^, ip, x), but will use rather another parametrisation formulated in 
spherical coordinates. This choice is made because it enables the efficient imposition 
of the boundary conditions in the numerical computations to follow. 

First we define the su(2) matrices in terms of Pauli matrices a a and spatial 
coordiates 9, ip 

= sin 9 cos nip 0\ + sin 9 sin nip a 2 + cos 9 a 3 , (24) 

Tg™' > = cos 9 cos nip o\ + cos 9 sin nip o 2 — sin 9 a 3 , (25) 

= — sin nip o x + cos nip a 2 , (26) 

where n is the winding number. For the Ansatze of the gauge field and the Higgs 
field we impose axial symmetry and in addition parity reflection symmetry. Using 
spherical coordinates the Ansatz of the gauge field 

Aidx 1 = A r dr + A e d9 + A^dip (27) 

can be parametrized by four functions Hk(r, 9), 

* = (1 _* M)) !£, (28) 

Ay = -nsin9 (H 3 (r,9) -f- + (1-H 4 (r,9)) -^-) . (29) 
The Ansatz for the Higgs field can be parametrised by two functions <j>k{r, 9) 

1 ?VH ?Y (n) 

-<P = Mr,0)-^ + Mr,0)^- (30) 



Expanding the curvature tensor in terms of the matrices ir^ /2, 

Fa = F{?^-, (31) 



we find for the non- vanishing coefficients F^\ 

4v) 1 



F%> = --(rdrHi + deHj , (32) 
n 



F^ = --sin 9{rd r H 3 - H ± H 4 ) , (33) 
= - sin 9{rd r H A + H x (H 3 + cot &)) , (34) 



r 



FtJ = -n sin 9{d e H 3 + H 3 cot 9 + H 2 H A - 1) , (35) 
F$ = n sin 9{8 e H A + (H 4 - H 2 ) cot 9 - H 2 H 3 ) . (36) 
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Similarly we find for the expansion of the covariant derivative of the Higgs field, 



the non-vanishing coefficients 

(A?0) (r) ■ 



r 

-(rd r (f)2 -Hxfa) 
r 



de<pi - H 2 4>2 , 
de4> 2 + H 2 0i , 

n sin 9 (fa H A + 2 (H 3 + cot 0) ) 



(37) 

(38) 

(39) 

(40) 
(41) 
(42) 



With these Ansatze the Hamiltonian densities eqs. (f|)-© can be expressed as 



+f£\d 9 <i>)U + fW{d 9 <i>)W 



sin 2 Or 4 ' 



(43) 



6 4 



sin 2 Or 2 
1 



2(^ - J(0? + 2 2 ))^ + (ZV) (9) (^0) ( ^ 2 



sin 2 #r 2 
1 



+ 
+ 



sin 2 9r 4 
1 

sin 2 6r A 



2(r/ 2 -^(0? + ^))F« + (^0)W(^0)^ 



(44) 



(2) 



54(rf 



2 + 2 )) 2 {[(^0)Ml 2 +[(D r 0)W 



(r) 



i 2 



(2) 



108 



+ 



+ - 



sin 2 Or 2 



(r) 



(45) 
(46) 



respectively. 
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With the Ansatz of the gauge fields we express the gauge constraint eq. (|2"3|) in 
the form 

G c = i(0 r i?! - ~d e H 2 ) = . (47) 

The variational equations are to be derived from the Hamiltonians of the models 
with the term [IlG 2 c added, where /i^ is the Lagrange multiplier. 

For n = 1 the spherically symmetric Ansatz is recovered with Hi(r, 8) = H 3 (r, 9) = 
4>2( r ,@) = and H 2 (r,9) = H^r, 6) = f(r) , <fii(r, 6) = h(r), where f(r), h{r) are 
defined in ||T5[ 



The boundary conditions follow from the requirements of finite energy and ana- 
lyticity as well as symmetry. They are given by 



at the origin 

#i(O,0) = # 3 (o,0) = o 

H 2 (0,9) = H 4 (0,9) = 1 

0i(O,0) = fa(P,0) = 



on the z-axis 

HfaO) = H 3 (r,0) = 

d e H 2 {r,0) = d e H 4 {r,0) = 

2 (r,O) = d e Mr,0) = 



at infinity 

H k (oo,6) = 

>i(oo,0)/27/ = 



on the p-axis 
#!(r,f) = H 3[ r 
d e H 2 {r,%) = d e H A 
i 2 (r,|) = 



' 2/ 



r > 2 J 



(4J 









4 Numerical results 



This Section deals with the results of the numerical integrations. The topological 
charge- 1 and charge-2 solutions of the systems under consideration will be found nu- 
merically and the energies per unit topological charge of these solitons will be calcu- 
lated as functions of the dimensionless coupling constants with the aim of deducing 
the presence of attractive and repulsive phases as the case may be, following the pro- 
cedure of Ref. [10]. This will be done by holding all but one of the dimensionless 



coupling constants constant and varying the last one only, although, given that our 
models are charcterised with more than one such constant, it would be possible to 
vary more than one of these at a time. This we do not do here. The present Section 
consists of three subsections. 

The first subsection is concerned with certain extended p = 1 models. Of these 
the first is the GG model, while the next five are perturbative extensions of the p = 1 
model which were described in Section 2.3.2 as perturbed p = 1 composite models 
consisting of the p = 1 model and all or some of the terms in the p = 2 model. The 
latter therefore differ from the GG model in two respects. Firstly, they incorporate 
terms from the p = 2 model which the GG model does not, and secondly, these 
terms may be regarded as perturbations of the p — 1 sytem unlike the GG model as 
explained in Section 2.3.2 and in Ref. 0. In particular, we have studied the p — 1 
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model extended by each of the four terms constituting the p = 2 model, separately. 
The interest in these (perturbative) models is that according to our analysis in section 
2.3.2 the attraction or repulsion of two like monopoles is decided by the slope of the 
energy versus the A a in question, and assuming quite reasonably that these curves are 
monotonic, this yields systems supporting only an attractive phase, or only a repulsive 
one. In addition to these four examples, we have studied the case of two nonvanishing 
\a, each chosen such that in the previous analysis the model with that particular A^ 
supported, respectively, an attractive and a repulsive phase. In such a case, it is 
intuitively expected that the combined terms will result in a model supporting both 
phases, and since both A^ are very small, the crossover point of the two phases is 
expected to be very near Xa = 0, the critical value of the model, which in this case 
coincides with the critical value of the PS model. 

The second subsection deals with the generic case of the composite model, i.e. 
that in which the values of the coupling constants A^ are not small, so that this 
is not a perturbative extension of the p = 1 model. We will set three of the four 
dimensionless coupling constants A^ equal to 1 and varying the fourth one. The main 
interest in these models is the question as to what phases are supported, the models 
themselves not being as closely almost self-dual as the p = 2 model on its own. For 
that reason it is not so easy to predict where, if anywhere, the crossover points will be 
situated. The other, technical, interest in the full composite model is that it affords 
us with an approach to the numerical integration of the p = 2 model on its own, by 
gradualy decreasing its p — 1 content. 

The third subsection deals exclusively with the p = 2 model, which from the 
viewpoint of numerical integrations is the most difficult of the models considered 
here. 

The partial differential equations are solved numerically using the iterative Newton- 
Raphson method p2~| . The calculations are performed with the rescaled Higgs field 
functions <f)i(r,8)/r], 02 (r, 6)/rj. We employed the compactified dimensionless radial 
coordinate x = rjr/(l + rjr) to map spatial infinity to the finite value x — 1. The 
differential equations are discretized on a non-equidistant grid in x and 9 covering 
the integration region < x < 1, < 9 < 7r/2. Typical grids have sizes of 70 x 20 
gridpoints. The error estimates are of the order 10~ 3 . 

4.1 Extended p = 1 models 
4.1.1 The Georgi-Glashow model 

This is the best known extended p = 1 model, and has been extensively studied in 
the literature analytically. A thorough numerical study of this model was recently 
performed in where it was confirmed that this system supports a repulsive phase 
only. We refer to || for the analysis of this model. 

This system consists of the p = 1 model to which the usual symmetry breaking 
Higgs self-interaction potential is added. It is well known || [J that the addition of 
the latter term cannot be regarded as a perturbation of the p = 1 model. In this 
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respect this model differs from the composite models following in the next sections. 

In Fig. 1 we show the energy density of the charge-2 solution of the p = 1 model 
as a function of the coordinates p = rjr sin 9 and z = rjr cos 9. The energy density 
has its maximum on the p-axis and is monotonic decreasing on the z-axis. Thus its 
shape deviates considerably form spherical symmetry. The surfaces of equal energy 
density form spheres for large distances form the origin and doughnuts for distances 
more close to the origin. The Maximum of the energy density corresonds is located 
on a circle in the x — y plane. 



4.1.2 A 4 = A 3 = A 2 = perturbed composite model 

In this model only the potential term 18x32 77~ 4 ^i'^i 2 ' > is added to the p = 1 model. 
In Fig. 2 the quantity AE(Xi), as defined in eq. ([TID , is plotted as a function of 
the coupling constant Ai (solid line a) for < Ai < 1. AE(Xi) is a monotonically 
increasing function. Consequently there is only a repulsive phase in this model. 

The dashed straight line in Fig. 2 corresponds to the linear approximation of AE 
defined in eq. (|20|) . The deviation of the function AE(Xi) from this line occurs at 
rather small values of Ai. This shows that the linear approximation is valid only for 
very small values of Ai, roughly in the range < Ai < 1/100. 

Note, that for this model the analysis of ref. H can be adopted. The reason is, 
that the actual form of the Higgs potential does not enter the analysis of the zero 
modes. Thus the conclusions of ref. M are valid for this model, too, provided the 
Higgs potential forces the Higgs field to decay faster than 1/r. 



4.1.3 A 4 = A 3 = Ai = perturbed composite model 

This model consists of the p = 1 model with only the term 18 x 32 r 7~ 4 ^27"4 2 ' > added to 
the p = 1 model. The quantity AE(X 2 ) is shown in Fig. 2 as a function of A 2 (solid 
line b). Again, AE(X2) is a monotonically increasing function. Thus there is only a 
repulsive phase in this model. 

The dashed straight line shows the linear approximation of AE, eq. fl20|) . For A 2 > 
0.1, AE deviates from this line. Therefore, the validity of the linear approximation 
to this model is restricted to the range < A 2 < 0.1 



4.1.4 A 4 = A 2 = Ai = perturbed composite model 

This model consists of the p = 1 model with only the term 18x32 ? 7~ 4 ^37'4 2 ' > added 
to the p — 1 model. This Skyrme like term leads to an attractive phase as can be 
seen form Fig. 2 (solid line c). AE(X^) is a monotonically decreasing function. The 
deviation from the linear approximation (dashed line) occurs at A3 = 0.1 implying 
that the validity of the linear approximation to this model is restricted to the range 
< A 3 < 0.1. 
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4.1.5 A3 = A2 = Ai = perturbed composite model 

This model consists of the p — 1 model with only the term j^^V" 4 ^^^ added. 
Again this model has only an attractive phase. The quantity AE(X4) (solid line d) 
as a function of A4 is shown in Fig. 2 together with the linear approximation (dashed 
line). The linear approximation is valid for < A4 < 0.2. 



4.1.6 A 4 = A 2 = composite model 

The numerical results for the composite models with only one of the terms eq. 
added to the p = 1 model show that there exist models with attractive or repulsive 
phases only. On the other hand, from the analysis of section 2.3.2 follows that for 
small coupling constants Xa models with both attractive and repulsive phases can 
exist, provided at least one of the terms added to the p = 1 model leads to a repulsive 
phase and the other one leads to an attractive phase. 

We studied numerically a model where the terms lg ^ 32 ?7~ 4 Ai?i^ 2 ' > and 18 x 32 ? 7~ 4 ^37"4 2 ' ) 
are added to the p = 1 model. In Fig. 3a we show the energy per unit topological 
charge E(Xi, X 3 )/n for the charge-1 and charge-2 solutions as functions of A3 with 
Ai = 0.01 fixed. Both functions increase with A3. For small values of A3 the energy 
per unit topological charge of the charge-2 solution (dashed-dotted line) is larger then 
the energy per unit topological charge of the charge-1 solution (solid line), whereas 
for larger values of A3 the energy per unit topological charge of the charge-1 solution 
is larger than the energy per unit topological charge of the charge-2 solution. The 
crossover occurs at A3 ~ 0.22 where both solutions have the same energy per unit 
topological charge E/n ~ 1.085 x Ann. Notice that attraction/repulsion occur for 
values of A3 above/below the crossover point, in contrast with the situation [10] in 
the Abelian Higgs model. 

As expected, this model supports an attractive and an repulsive phase. The 
crossover is at a small value of the coupling constants Xa and the energy per unit 
topological charge at this point is near the Bogolmol'nyi bound. However, in the ex- 
ample chosen the crossover is outside the range where the linear approximation is 
valid. From the linear approximation we can only predict the existence of attrac- 
tive/repulsive phases, but not the features of the models in detail. 

Another example is shown in Fig. 3b. Here the energy per unit topological charge 
is shown for the charge-1 solution (solid line) and the charge-2 solution (dashed- 
dotted line) as a function of Ai for fixed coupling constant A3 = 1. Again, for both 
solutions the energy per unit topological charge is an increasing function. There is an 
attractive phase for small values of Ai and a repulsive phase for large values of Ai. 
The crossover is at Ai ~ 0.07 where the energy per unit topological charge of both 
solutions is E/n = 1.375 x 4nn. For large values of the coupling constant, like A3 = 1, 
the linear approximation is certainly not valid anymore. Still the model supports 
both an attractive and a repulsive phase. However, the energy per unit topological 
charge at the crossover now exceeds the Bogolmol'nyi bound considerably. 
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4.2 Composite model 

In this section we consider the model defined by the Hamiltonian eq. (0), where the 

parts of the p = 1 model and the p = 2 model are treated on the same footing. The 

linear approximation of the energy per unit topological charge as a function of the 

coupling parameters A^ is not applicable for this model. 

As a typical example we discuss the charge-2 solution for the coupling constants 
= 1, A — 1, . . . , 4. The energies per unit topological charge of the charge-2 and 

charge-1 solutions are E n= 2/2 = 2.1433 x Attt] and E n= \ = 2.05993 x Airrj, respectively. 

Thus the charge-2 solution is in the repulsive phase for these values of the coupling 

constants. 

In Figs. 4a-4g we show the energy density and the field configuration for the 
charge-2 solution as functions of the dimensionless coordinates p and z. In Fig. 4a the 
energy density is shown. It possesses a maximum on the p-axis and it is monotonically 
decreasing along the z-axis, indicating a strong deviation form spherical symmetry. 
Like for the p — 1 model the surfaces of equal energy density form spheres und 
doughnuts. 

Let us now compare the energy densities of the charge-2 solutions in the composite 
model and in the p = 1 model, shown in Fig. 1. At the origin the energy density 
of the composite model is roughly one order of magnitude larger than the energy 
density of the p — 1 model. In the composite model the energy density approaches 
its asymptotic value at a shorter distance form the origin than in the p = 1 model. 
Thus the composite model possesses a core with a high energy density, whereas the 
core in the p = 1 model is softer and more extended. 

In Fig. 4b the gauge field function H\ (p, z) is shown. According to the boundary 
conditions Hi(p, z) vanishes on the p and z axis. It has a maximum located near the 
origin. For spherically symmetric solutions Hi(p, z) would vanish identically. Thus 
Hi (p, z) is excited due to the axial symmetry of the solution. Fig. 4c shows the gauge 
field function H2(p,z). The deviation from spherical symmetry for this function is 
small near the origin and nearly vanishing far away from the origin. The function 
H 3 (p, z) is exhibited in Fig. 4d. Simmilarly to the function Hx(p, z) shown in Fig. 4b 
this function is excited due to the axial symmetry of the solution. Its maximum is 
located near the origin. Note that it is roughly one order of magnitude smaller than 
the maximum of the function Hi(p,z). In Fig. 4e the function Hi(p,z) is shown. 
Like the function H2(p,z) shown in Fig. 4.c the deviation from sphercial symmetry 
is only small for this function. The Higgs field function <fii(p,z) is shown in Fig. 4f. 
At the origin this function is linear in z for p = whereas for z = it is quadratic in 
p. For large distances from the origin the deviation from spherical symmetry almost 
vanishes. Fig. 4g shows the Higgs field function <p2{p,z). It has a minimum near 
the origin and vanishes on the p and z axis due to the boundary conditions. Like 
the gauge field functions Hi(p,z) and H 3 (p, z) this function is excited by the axial 
symmetry of the solution. 

The gauge field functions and Higgs field functions possess a non-trivial depen- 
dence on the coordinates p and z only in the region near the origin where the core 



15 



of the energy density is located. As a result of our numerical integrations we have 
found that outside this core the functions decay to their vacuum values at infinity at 
a rate faster than a power behaviour. This is in contrast to the solutions of the p — 1 
model, where in absense of the Higgs potential the Higgs fields possess a power law 
decay and only the gauge field functions decay exponentially. Thus for the composite 
model the asymptotic decay of the energy density outside the core originates not from 
the Higgs field functions but from the gauge field functions only. 

Let us next look at the dependence of the energy on the coupling constants A^ 
for the charge-1 and charge-2 solutions of the composite model. 

We have solved numerically the partial differntial equations for fixed A4 = 1 and 
several values of A^ in the range < A^ < 1.5 for A = 1,2,3. The results are exhibited 
in Fig. 5. Lines a and a' represent the energy per unit topological charge of the charge- 
1 and charge-2 solutions, respectively, as functions of the coupling constant Ai with 
A2 — A3 = 1 fixed. The lines cross at Ai ~ 0.1. Thus there is both an attractive 
phase (Ai < 0.1) and a repulsive phase (Ai > 0.1). The value of the energy per unit 
topological charge E/n = 1.865 x Attt] at the crossover is rather large compared to the 
Bolgolmol'nyi bound, which takes the value 1.1 x Atti] for Ai = 0.1 in this case. 

Lines b and b' represent the energy per unit topological charge of the charge-1 
and charge-2 solutions, respectivley, when A2 is varied with fixed Ai = A3 = 1. In the 
range < A2 < 1.5 the lines b and b' are almost parallel and there is no crossover. 
The upper lines correspond to the charge-2 solutions. There is a repulsive phase only 
in this case. 

Finally we fix Ai = A2 = 1 and vary A3. The energy per unit topological charge 
is represented in Fig. 5 by lines c and c' for the charge-1 and charge-2 solutions, 
respectively. Again, lines c and c' are almost parallel. Thus, there is only a repulsive 
phase. 

4.3 p = 2 model 

The numerical integration of the variational equations of the p = 2 model is a very 
difficult task. In the absense of the terms of the p = 1 model numerical instabilities 
arise leading to a failure of the algorithms. We succeded to construct the spheri- 
cally symmetric charge-1 solutions. However, we did not find the axially symmetric 
charge-2 solutions of the p = 2 model on its own. Instead we used an approxima- 
tion procedure to calculate the energy of the charge-2 solutions as explained in the 
following. 

We start with the composite model and regard now the p = 1 part as a perturba- 
tion. For that purpose we modify the static Hamiltonian eq. ( |i~3"l) by introducing the 
new coupling constant Ao, 

n {c) = x n {1) + r]^n {2 \ (49) 

and consider the limit Ao — > 0. Assuming the existence of the solutions to the p = 2 
model, we will now consider the model (^) as a perturbation of the latter by the 
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p = 1 model, characterised by a small value of the dimensionaless coupling constant 
A < 1. 

We expand the energy of the magnetic charge-n soliton in the vicinity of the 
parameter Ao, 



£ (n) (A + AA ) = £ (n) (A ) + AA ( 



dE^ 



d\ 



+ o(\ 2 ). (50) 

Ao 



dXo 



Ao=0 



can be calculated 



With the same reasoning as in section 2.3.2, the quantity 

exclusively in terms of the exact solutions. Now, because we assumed that Ao is very 
small, we can evaluate ([50]) at AAq = — Aq. Thus for sufficiently small Ao, 



£ (n) (0) « £ (n) (A ) - A ( 



d\ 



(51) 



will give a very good approximation to the energy of a solution of the p = 2 model. 
This energy consists of two terms, the first is the energy of the modified p = 2 model 
for a small value of Ao- The second term is given by — Ao / 7i^d 3 r evaluated for the 
solution with finite (but small) Ao- 

We can now estimate the energy of the p = 2 model on its own by calculating the 
energy of the modified composite model for a very small value of Ao and by subtracting 
the contribution of the p — 1 part, i. e. the quantity Ao / H.^d 3 r evaluated with the 
solution for the small value of Ao- 

To check the validity of our limiting scheme, we followed this procedure for a set 
of decreasing values of A and found, that the estimated energy becomes indeed a 
constant for Ao small enough. 

Let us first we discuss the properties of charge-2 solutions for the p = 2 model at 
the critical values of the coupling constants A^ = 1, A = 1, . . . , 4, determined by the 
investigation of the charge- 1 solutions. The energy per unit topological charge of the 
charge-1 solution was found to be E^ 1 ' (\a = 1) = 1.00 02 x 47r^, whereas for the charge- 
2 solutions we find for the same coupling constants E^ 2 \\a = l)/2 = 1.003 x Awq. 
Thus the like monopoles repel for the critical values of the coupling constants. If 
we compare with the Bogolmol'nyi bound, given by Anr], we see that for both the 
charge-1 and charge-2 solutions the energies per unit topological charge are very close 
to the lower bound, indeed. 

In order to compare the p = 2 model with the p = 1 model and the composite 
model we show in Fig. 6a-6g the energy density, the gauge field functions and the 
Higgs field functions for the solution of the modified model (f49|) for the coupling 
constants Ao = 0.01, Xa = 1, A = 1, . . . , 4 and assume that this configuration is close 
to the solution of the p = 2 model on its own. 

The energy density, shown in Fig. 7a as a function of the dimensionless coordinates 
p and z, is similar to the energy density of the composite model, shown in Fig. 4a. It 
possesses a maximum on the p-axis near the origin and it is monotonically decreasing 
along the z-axis. Its value at the origin is of the same magnitude as for the composite 
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model and the decay for large distances from the origin is slightly faster. Like for 
the p = 1 model and the composite model the surfaces of equal energy density form 
spheres und doughnuts. 

The gauge field functions and the Higgs field functions are shown in Figs. 6b-6e 
and Figs. 6f-6g, respectively. They look very simmilar to the solution of the composite 
model shown in Figs. 4b-g, except for the gauge field function H 3 (p, z). For the p = 2 
model H 3 (p,z) is negative and its minimum is of the order « —0.1, whereas for the 
composite model it is positive with a maximum of the order « 0.025. 

Let us next look at the dependence of the energy on the coupling constants A^ for 
the charge- 1 and charge-2 solutions. For the p = 2 model we have fixed the coupling 
constant A4 = 1. For the charge-1 solutions we solved the ordinary differential equa- 
tions of the p = 2 model for several values of the coupling constants. The energy for 
the charge-2 solutions are calculated approximately as discussed above. 

The results are shown in Fig. 7a-7c. In Fig. 7a we exhibit the energy per topological 
charge of the charge-1 solution (solid line) and the charge-2 solution (dashed-dotted 
line) as a function of the coupling constant Ai with fixed A2 = A3 = 1. Both lines cross 
at Ai = 0.89. There is an attractive phase for Ai < 0.89 and a repulsive phase for 
Xi > 0.89. The energy per topological charge at the crossing is E/n = 0.98 98 x Aitn. 
As expected the crossover takes place at values of the coupling constants near the 
critical values Ai = A2 = A3 = 1. 

In Fig. 7b we exhibit the energy per topological charge of the charge-1 solution 
(solid line) and the charge-2 solution (dashed-dotted line) as a function of the coupling 
constant A2 with fixed Ai = A3 = 1. In this case the crossing is at A2 = 1.07, with a 
repulsive phase for A 2 < 1.07 and an attractive phase for A 2 > 1.07. For the energy 
per topological charge at the crossing the value E/n = 1.0333 x 4^77 was found. Again, 
the values of the coupling constants at the crossing are near the critial values. 

In Fig. 7c we exhibit the energy per topological charge of the charge-1 solution 
(solid line) and the charge-2 solution (dashed-dotted line) as a function of the coupling 
constant A3 with fixed Ai = A2 = 1. In this case the crossing occurs for a value of 
the coupling constants A3 = 1.49, which is considerably larger then the critical value 
A3 = 1. There is a repulsive phase for A3 < 1.49 and an attractiv phase for A 3 > 1.49. 
The energy per topological charge at the crossing takes the value E/n = 1.1969 x 4 pin. 

Notice that the situation in Figs. 7b and 7c, namely that of attraction/repulsion 
taking place for values of the coupling constants larger /smaller than that of the 
crossover point, contrasts with the corresponding situation [110] for the Abelian Higgs 
model. 

5 Summary and Conclusions 

The main aim of the present paper was to acheive monopole field configurations in 
which like-monopoles can attract as well as repel. We found this to be the case in var- 
ious new 3 dimensional 50(3) gauged Higgs models with isovector Higgs fields. Such 
a property contrasts with the 't Hooft-Polyakov monopoles of the Georgi-Glashow 
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model, in which like-monopoles only repel. The main model we employed was one we 
described as the p = 2 model because it is the second one in a hierarchy of models 
descended from 4p dimensional YM systems ||13|| , the first one of which, p — 1, is the 
GG model in the PS limit. The spherically symmetric magnetic charge-1 solutions 
of the p = 2 model were previously studied, from which we concluded that the type 
of behaviour sought in the present work was likely to materialise. Accordingly in 
the present investigation, we studied the charge-2 axially symmetric solutions of this 
model, and found indeed that for some values of the dimensionless coupling param- 
eters the interaction energies of the monopoles led to bound states. In this sense, 
the 3 dimensional p = 2 model is very similar to the 2 dimensional Abelian Higgs 
model, both featuring a crossover between an attractive and a repulsive phase. In- 
deed, our procedures were the same employed by Jacobs and Rebbi JED] in the 2 



dimensional case. There is however, a qualitative difference between these models in 
the sense that the attractive and repulsive phases of the p = 2 model occur on the 
opposite sides of the crossover points for some of the coupling constants. Another 
difference between these models is that in the p = 2 case there is really no genuinely 
non-interacting, minimal, field configuration. Nevertheless there is an almost non- 



interacting field configuration which was discussed in some detail in | 15| , which turns 
out to be approximately the crossover point between the two phases. 

In addition to the p = 2 model, we studied combinations of the p = 1, and the 
p = 2 models, as well as extensions of the p = 1 model by some of the constituent 
terms in the p = 2 model which is tantamount to adding a particular Skyrme-like 
term. In some of the latter examples we also considered small perturbations of the 
p = 1 model, something which is radically different from adding the Higgs potential 
to the p = 1 model to obtain the GG model. ( Recall that the BPS solutions cannot 
be related to the monopoles of the GG model perturbatively.) Perhaps the most 
interesting outcome of the latter considerations was the construction of extended 
p = 1 models whose monopole solutions feature either an attractive phase only or a 
repulsive phase only, between two like-monopoles. 

Apart from their intrinsic interest in providing generalised monopoles exhibiting 
both attractive and repulsive phases, it is hoped that this property of the gener- 
alised monopoles may be of some potential interest in the role of monopoles in grand 



unified theory, and possibly to the cosmic monopole problem [23]. The models we 
have considered could occur as extensions of the usual models of the grand unified 
theory. In such models there is attraction between like-monopoles (and ofcourse like- 
antimonopoles) as happens with type I cosmic strings p4| . We do not speculate here 
on the possible consequences of this analogy. Perhaps the most likely effect of such 
(generalised) monopoles may occur in the gravitational clumping [^] of monopoles 
leading to the annihilation of monopole-antimonopole pairs. In contrast with the 
usual theory, there will be only attractive forces between the monopoles, indepen- 
dently of the gravitational interaction. 
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Figure captions 



Fig. 1 

The energy density in units of rf of the p = 1 model is shown as a function of the dimen- 
sionless coordinates p = rjr sin 6 and z = r/r cos 9. 

Fig. 2 

The quantity AE(Xa) = E < - 2 \Xa)/2 — E^\Xa) in units of Airr/ (solid lines) is shown as a 
function of with Xb = 0, B ^ A. The dash-dotted lines indicate the slope at Xa = 0. 
Line a corresponds to Ai / 0, line b to A2 / 0, line c to A3 / and line d to A4 / 0, 
respectively. 

Fig. 3a 

The energy per unit topological charge in units of Attt] is shown as a function of A3 with 
Ai = 0.01, A2 = A4 = for the charge n = 1 solution (solid line) and for the charge n = 2 
solution (dash-dotted line). 

Fig. 3b 

The same as Fig. 3a for Ai with A3 = 1 and A2 = A4 = 0. 
Fig. 4a 

The energy density in units of r\ of the charge n = 2 solution of the composite model with 
\a = 1, A = 1, . . . ,4 is shown as a function of the dimensionless coordinates p = i]rsm9 
and z = rjr cos 6. 

Fig. 4b 

The same as Fig. 4a for the gauge field function H\{p, z). 
Fig. 4c 

The same as Fig. 4a for the gauge field function H2(p, z). 
Fig. 4d 

The same as Fig. 4a for the gauge field function H^(p, z). 
Fig. 4e 

The same as Fig. 4a for the gauge field function H^p, z). 
Fig. 4f 

The same as Fig. 4a for the Higgs field function 4>i(p, z). 
Fig. 4g 

The same as Fig. 4a for the Higgs field function (p, z). 
Fig. 5 

The energy per unit topological charge in units of 4irr] of the composite model is shown as 
a function of Xa, A = 1,2,3 with fixed Xb = 1, B 7^ A for the charge n = 1 solutions 



22 



(solid lines) and for the charge n = 2 solutions (dash-dotted lines). Lines a and a' corre- 
spond to A2 = A3 = A4 = 1, lines b and b' to Ai = A3 = A4 = 1, and lines c and c' to 
Ai = A2 = A4 = 1, respectively. The thin solid line indicates the Bogolmol'nyi bound. 

Fig. 6a 

The energy density in units of rf of the charge n = 2 solution of the p = 2 model with 
A,4 = 1, A = 1, . . . ,4 is shown as a function of the dimensionless coordinates p = r\r sin 6 
and z = r/r cos 9. 

Fig. 6b 

The same as Fig. 6a for the gauge field function H\(p, z). 
Fig. 6c 

The same as Fig. 6a for the gauge field function H2(p, z). 
Fig. 6d 

The same as Fig. 6a for the gauge field function H^(p, z). 
Fig. 6e 

The same as Fig. 6a for the gauge field function H±(p, z). 
Fig. 6f 

The same as Fig. 6a for the Higgs field function 4>i(p, z). 
Fig. 6g 

The same as Fig. 6a for the Higgs field function foip, z). 
Fig. 7a 

The energy per unit topological charge in units of Attt] of the p = 2 model is shown as a 
function of Ai with fixed = 1, A / 1 for the charge n = 1 solution (solid line) and for 
the charge n = 2 solution (dash-dotted Line). The thin solid line indicate the topological 
bound. 

Fig. 7b 

The same as Fig. 7a for A2 with fixed A^ = 1, A ^ 2. 
Fig. 7c 

The same as Fig. 7a for A3 with fixed A^ = 1, A ^ 3. 
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1/n dE/dX A \x A=0 


CA 


A 


n = 1 n = 2 




1 


0.688 2.925 


2.237 


2 


0.383 0.437 


0.054 


3 


0.367 0.291 


-0.076 


4 


0.061 0.020 


-0.041 



Table 1 

The quantity 



1 dEW 



n d\A 



X A =0 



defined in (|T^) for the Hamitonians Ti^ is given for the charge- 1 



and charge-2 solutions together with their difference ca in units of 4nrj. 
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Fig. 1 
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Cl 

to* 
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Fig. 3a 

1 .12 1 1 
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Fig. 3b 

1.6 i i 
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Fig. 4a 
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Fig. 4c 
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Fig. 4e 
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Fig. 4f 
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Fig. 4g 
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Fig. 6a 
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Fig. 6c 
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Fig. 6e 
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Fig. 6f 
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Fig. 6g 
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Fig. 7a 



1.005 



0.995 



LU 



0.99 



0.985 



0.98 



0.975 
0.8 
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Fig. 7c 
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